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Abstrac t 

Tire development of control of large flexible structures 
technology must include practical demonstrations to aid in 
the understanding and characterization of controlled struc- 
tures in space. To support this effort, a testbed facility 
has been developed to study practical implementation of 
new control technologies under realistic conditions. The 
paper discusses the design of a second order, acceleration 
feedback controller which acts as an active vibration ab- 
sorber. This controller provides guaranteed stability mar- 
gins for collocated scnsor/actualor pairs in the absence of 
scnsor/actuafor dynamics and computational time delay. 
Experimental results in the presence of these factors are 
presented and discussed. The primary performance objec- 
tive considered is damping augmentation of the first nine 
structure I modes. Comparison of experimental and pre- 
dicted closed-loop damping is presented, including test and 
simulated time histories for open and closed-loop cases. 
Although the simulation and test results arc not in full 
agreement, robustness of this design under model uncer- 
tainty is demonstrated. Ibc basic advantage of this second- 
order controller design is dial the stability of the controller 
is model independent. 

Introduction 

l arge space structure control requires special control 
design methodology. High modal density in the controller 
bandwidth makes rolloff of control authority a problem. 
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Many of the theoretical issues with regard to large space 
structure control were recognized and initially treated by 
Benhabid, Ftashner, and Tung . 1 Experimental validation is 
now being demonstrated with a series of facilities which 
have been designed for that purpose . 2 ' 5 ' 4 Ft is recognized 
that the complexity of these structures foster the use of con- 
trollers with high and low authority. Theoretically, low au- 
thority controllers arc used to augment the inherent damp- 
ing of the structure, while the high authority controllers arc 
designed to achieve the desired performance. One form of 
these low authority controllers is an active vibration ab- 
sorber (AVA), which is the subject of this paper. 

NASA and the Controls Structures Interaction (CSI) 
program have undertaken the design, development, and 
fabrication of generic structures to study fundamental prob- 
lems in the Implementation of advanced control method- 
ologies. The program stresses the application of existing 
sensor and actuator technologies for identification and con- 
trol. The first such testbed of the CSI program has been 
completed at the NASA Langley Research Center. The 
configuration selected, which is described in detail in a later 
section, uses thrusters and accelerometers as its primary 
actuation and sensing devices. Initially, controller designs 
that would not require knowledge of the dynamics of the 
structure are preferred. This allows the implementation of 
such controllers even when model uncertainties arc large. 
For this purpose, a procedure for designing second order 
controllers using passitivity concepts is discussed for con- 
tinuous time systems . 6 Using this approach, decentralized 
controllers have been designed and used for experimental 
validation. 

The paper outline is as follows: first the controller de- 
sign theory is presented followed by a mechanical analogy. 

A brief description of the testbed is given, and a controller 
design procedure for the laboratory model is described. Fi- 
nally, lest results from the experimental implementation are 
provided and compared with the analytical results. 
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Controller Deflign Thnnry 

For large scale systems, a model of the structure is 
orten obtained via finite element methods resulting In a set 
of state and output equations of the form: 

Mz »- Di f Kz — Du (|) 

y ~ fin S -f Il v i + U^x (2) 

Mere z is an n x i state vector, and M, D, and K are 
mass, damping, and stiffness matrices, respectively, which 
generally arc symmetric and sparse. The nxp influence 
matrix B describes the actuator force distributions for the 
^ x 1 in r ,,t vcc lor u. Equation (2) is a measurement 
equation with y as the m x l measurement vector, H a the 
m x n acceleration influence matrix, //„ the m x n velocity 
influence matrix, and H d the m x n displacement influence 
matrix. F< r direct output feedback control, the input vcclor 
« can he v/ritten 

1 ~ Gy = —GHaX — GH v i — GHjz (3) 

where G is a gain matrix to he determined. Substituting 
Vi}. (3) into Eq. (I) yields 

(M + BGH a )x + (D f BGH„)i + ( K f BGH d ) z = o (4) 


Let the input vector « in Fq. (I) ami in Fq. (5) be 
* — Ve = Hade + Hvdc + H df te (7) 

»< - » = Hat + H v i + H d z (g) 

Substituting Fq. (7) into Fq. (I) nnd Fq. (R) lnfr> Fq. (5) 
yields 

Af,5, t n t i t f K t z t = o (9) 

where 



If the design parameters arc chosen such that M t , D t and 
Kt are positive definite, the closed loop system, Eq. (9), 
becomes asymptotically stable. 

Consider the special ease for acceleration feedback only 
where H d — //„ — H de — ff vc — 0 in Eq. (9). This makes 
the closed-loop mass, damping and stiffness matrices 


It is shown by Junng and Phan 6 that for a structural sys- 
tem, output velocity feedback with collocated sensors and 
actuators makes the closed-loop system asymptotically sta- 
ble with an infinite gain margin. Without velocity sensors, 
however, the system damping cannot be augmented with 
ilircct output feedback unless additional dynamics arc In- 
troduced. 

Let the controller to be designed have a set of second- 
order dynamic equations and measurement equations sim- 
ilar to the system equations, Eqs. (I) and (2), 



For M t to be a positive definite matrix, it must be a real, 
symmetric matrix satisfying 

( 10 ) 


Md c f D e i e + K e x e =- B e u e (5) 

Ve ~ H a< .ic + H„de + H dc x e (6) 

Here x c is the controller state vcclor of dimension n e , 
and M e , f) e , and K e can be interpreted as the controller 
mass, damping, and stiffness matrices, respectively, which 
in general, arc symmetric and positive definite to make the 
controller asymptotically stable. The n e x m influence ma- 
trix He describes lire force distributions for the my 1 input 
force vector u e . Fquation (6) is the controller measurement 
equation having y e as the measurement vector of length p, 
H ae the fxn, acceleration influence matrix. Hue the px» e 
velocity influence matrix, and the p x n e displacement 
influence matrix. All the quantities u e , y e , and n e , are 
arbitrary which means that M e , D e , K e , H de , H„ t , Hae 
and B e arc the design parameters for the controller. 


for any real vector *, except the null vcclor lb make M, 
symmetric, it is required that 

BHae- hJbJ ( 11 ) 

Substituting the definition of z, and M t into Eq. (10) and 
using Eq. (1 1 ) yields 

= * T (M - BHaeflJ c B T )x+ 

{HI'B T z - Ze) T (Hj t B T x - Xe) + xJ(Me - I)xe 

This equation is greater than zero if H ae and M e arc chosen 
such that and M c — I arc positive definite. 

Note that this is a sufficient condition but not a necessary 
one. 
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In an attempt to obtain a controller structure without 
this restriction oil H ac and M e , let the u in Eq. (7) be 
modified to include a direct acceleration feedback term. 
Then 


which makes 


I >e lining 


|fc - Call - Hoc Xe — Gay 

(12) 

\M t nG a H a -niloc 1 

[ Delia Me J 

03) 

G a = HaeM~ l D e 

(14) 


Me,D e and K e to be positive definite and that Eq. (11) is 
satisfied. 


Mechanical A nal ogy 

With further constraints on the plant and controller, 
this second order controller can be thought of as attaching 
an additional mass, spring and damper assembly to the 
structure at the location of the sensors/actuators. For the 
ease of a single mode plant, a diagram representing the 
closed loop system can be drawn as shown in Fig. 3. 
Letting x e = x a - x results in plant and controller equations: 


substituting into Fq. (13). and using the symmetry condi- 
tion in Eq. (II). the positive definite condition in Oq. (10) 
Ik'coiiics 

if MfZt = x^Mx- 1- 

(Mc l n c ll a z - x c ) r M c (Mi l HclUx - x e ) (15) 

which is positive if M and M c arc positive definite ma- 
lices. Figure ! shows a block diagram of llie closed- 
imp system with acceleration feedback. Figure 2 depicts 
a reduced block diagram for the same system. To de- 
sign a second-order controller having the structure given 
in (5),(6),(8),(12), and (14), it is only necessary to require 



I • gure I. »Iock diagram of closed loop system with 
direct, and dynamic feedback 


Plant : Mx f Dx \ Kx - 

V - s * (16) 

Controller : M c z c t D c z c t K c r c ~ Mey 

Ze ~ K c z c DeZe (17) 

A block diagram for this physical system is shown in Fig. 
4. Upon comparison of Fig. 2 and Fig. 4, H ac and B c 
must satisfy 

little Me 1 (18) 



Figure 3. Single mode struct ure-rout roller design 
model 



figure 2 lleduced block diagram of closed loop controller closed loop system 
system with general acceleration feedback 











Assuming collocated sensors and actuators ami unity force the following section, a description of the testbed fac ility 

input, the constraints given in Fqs. (I I) ami (18) dictate: is presented. 1 


1 (19) 

n - 1 (20) 

//(ic = Me (21) 

tie = Me (22) 


Using the controller equations given in Rq. (17), a state 
space model can be written: 


i — Ax I- By 
ti - Cx f Dy 

where 


(23) 

(24) 


A - f ° 1 

[ Kef Me Pef Me 

C=\-K e -De I, D — (0) , S-r, |**J 

Ihc controller low frequency gain equals M e . In the pres- 
ence of sensor bias, a large M c could result in large ac- 
tuator command offsets or even saturation. It is important 
to select conservative values for Me to reduce lire cITccts 
of offsets. In analogy with vibration absorbers 0 , Ihc fre- 
quency of the absorlvcr is set equal to the targeted mode, 
i.c. u 2 - KefMe . After selecting values for M e and K e , 
values for D c arc selected based on root locus design. In 




LaboratoryModfiLDcscriirtioji 

A photograph of lire evolutionary model testbed (4J 
is shown in Fig. 5. Ihc structure weighs 741 lbs. and 
is supported from the ceiling with two sled cables with 
lengths of 64.5 ft. each. The truss is made of aluminum 
struts forming 10 in. cube bays of single lace diagonals 
alternating so that opposite faces cross. The batten plane 
has diagonals which also alternate. The major components 
of the structure arc the center section 52.5 ft. long, a 
16 ft. diameter reflector, and a 9.2 ft. tower where a 
laser beam is located. Ihc steel cables arc attached to the 
ends of two cross mendrer trusses 16.7 ft. long to ensure 
stability when suspended. Fire actuators are proportional 
thrusters using a 125 psi external air supply connected to 
the model via flexible rubber hoses. A total of R actuator 
pairs arc placed at four locations along the truss as shown 
in Fig. 6, and each actuator pair is capable of producing 
4.4 lbs. In addition, 8 accelerometers arc located with 
the thrusters for identification and control experiments. 
The finite element model has Rl modes below 50 lb. 
Because of the suspension, the first six structural modes 
arc pendulum modes with frequencies between 0.1 lb. and 
0.9 Hz. The next three modes are bending modes with 
frequencies between 1.4 Hz and 1.9 Hz. Orthogonal views 
of Ihc mode shapes Tor modes 6 through 9 arc provided 
in Figs. 7-10. Control laws arc implemented on a digital 
computer at sampling rales of 80 Hz ami 350 Hz. 



Reflector 


Figure 0. Srlirmatir <>f Hm f'SI Hvolnt ionary Model 
snowing nrhmlnr/qf tisf»r local im»q 
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Controller Design 

The design theory discussed is now used to develop 
.SISO controllers for each actualor/scnsor pair. The design 
approach is to determine the most control lablc/obscrvablc 
modc(s) for a particular input/output pair ami tune a sec- 
ond mder controller In die modc(s). Radi second order 
controller has three design parameters, namely, M e ,D e , 
and K t . If one actualor/scnsor pair is to control n sepa- 
rated modes, n distinct sets of design parameters need to 
lie determined. The procedure used to design controllers 
Tor the evolutionary model is the following: 

(1) Choose M e to avoid large command offsets. 

(2) Choose a Large! mode for each actualor/scnsor pair 
ami determine the value of K e needed to match its 
frequency, where K e /M e - w 2 . 

( d) Determine the value of De which optimizes the damp- 
ing value of the target mode. 


7 lie design objective is to increase damping in the first 
nine modes. A value of 0.1 is selected for M e based on 
previous experimental rcsulLs. The modes selected and the 
corresponding actualor/scnsor pair used to control tlicm 
are presenter! in table I. Each acluator/sensor location 
is labeled with a number from I to R, (sec Pig. f>). Using 
fable I, actualor/scnsor pair at location 4 is responsible 
lor controlling mode number 8. Tltc actuators at locations 
1.2,3, and 7 have two independent controllers for the 1** 
and 2 target modes, while the remaining four have one. 


Although the acluator/sensor at fix ation 8 is responsible for 
modes 6 and 7, only one controller is designed by tuning 
its frequency in between the two modes. To optimize 
damping, a root locus is performed for each actualor/scnsor 
pair by fixing M c and K e and varying D e . In locations 
where more than one controller is used, each controller 
is optimized separately. Figure 1 1 shows the r<x>t fix: us 
for mode 6 using actualor/scnsor pair 8 for values of D e 
ranging from 0.1 to 1.0. A maximum damping value of 
3.5% is achieved when /7 c -0.4 The controller parameter 
values, Mc,D f , and K e , for each second order controller 
designed are given in 75iblc 2. 
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Table I: Mode number to be controlled and corresponding acluator/sensor location 


Location No. 

l* 1 Target Mode 

2 n< * Target Mcxlc 

I 

3 

9 

2 

4 

8 

3 

2 

9 

4 

8 


5 

I 


6 

3 


7 

1 

8 

8 

6 

7 


Table 2: Controller design parameters 



l* 1 Ilirget Mode 

2 nJ Thrget Mode 1 

Location 

Me 

D e 

K c 

Me 

D e 

Ke 

1 

- 1 

.06 

.089 

.1 

.4 

1 3.22 

2 

.1 

.2 

2.10 

.1 

.5 

10.88 

3 

.1 

.04 

.083 

.1 

,4 

13.22 

4 

.1 

.4 

10.88 

, 


_ 

5 

.1 

.03 

.083 

- 


_ 

6 

.1 

.04 

.083 

- 

_ 

_ 

7 

.1 

.03 

.083 

J 

.4 

10.88 

8 


.4 

3.95 

- 

- 

- 
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Experimental Itesultg 

1. Implementation Iasuep 

TIic controllers described in the previous section have 
special properties lliat must be considered in the actual 
implementation. The active implementation of the AVA 
controller should ideally emulate a physical spring-mass- 
damjicr attached to a stmeture at a single point. This state- 
ment implies die critical assumption dial collocated sensors 
and actuators can be realized in the active implementation 
of the controller. Sensor and actuator collocation must 
he spaiial and temporal. Spatial collocation is difficult to 
achieve Ixxausc of ihc finite dimensions of sensors and ac- 
tuators. "typical mounting of sensors and actuators result 
in spatially close arrangements, not collocated. Temporal 
collocation is even more difficult to achieve. First, the sen- 
sors and actuators must have no inherent dynamics in the 
frequency range zero to infinity. Second, the control law 
calculation must produce no phase shift except for that pro- 
duced by the AVA controller, lienee, temporal collocation 
can be most closely achieved with very large bandwidth 
sensors and actuators and analog computations. 

In Hie following paragraphs, die digital implementation 
of the AVA controller will tie demonstrated on the structure 
shown in Fig. 5. Tltc effects of sample delay, the sequence 
of real time computations, and spatial non collocation arc 
shown using experimental data. It should be noted that the 
sensors used in this study have a bandwidth of 300 Ilz, 
whereas the actuators have a bandwidth of about 45 Hz. 
In our ease, the maximum frequency to be controlled is at 
1.9 II/, and thus, the assumption of negligible sensor and 
actuator dynamics is not grossly violated. 

Figure 12 shows the closed-loop implementation of the 
system controller. 7Tic external disturbance r produces the 
structural acceleration output y. Switches A and D arc 
used to close the controller loop. Classic implementation 
of Kalman filter based controllers call for switch A to 
l>c closer! at all limes since the controller is actually an 
estimator of die structural slates. It will be shown that 
die AVA controller is not an estimator and switches A 
and B should be dosed simultaneously lo activate closed- 
loop control. When closed-loop control is activated the 
following steps bike place at each of die k ih sample 
perirxls: 

( 1) Input y(k) 

( 2 ) Out|HJ( u(k) 

(3) Compute control force u(it \ 1) = Ci(k) + Dy(k) 

(4) Compute controller slates *(* * l) = Qz(k) + ty~(k) 

(5) Wait for next clock pulse 


where and V are Ihc time discrcli/cd versions of die A 
and B man ices in Eq. (23). It is also noted that die D 
matrix for die AVA controller is null. 


The AVA controller has been implemented following 
the sequence of steps above at a sample rale of 80 Hz. 
The structure was excited at two bending modes and 
two pendulum modes for 10 seconds, then closed-loop 
control was initiated at 12.5 seconds. Figure 13 shows the 
acceleration output with switch A closed from 1=0 seconds 
and implementing steps 1-5 above in sequence. 7T»c 
controller is unstable in a mode near 7 Hz. Since die AVA 
controller is not an estimator, switch A and B were closed 
simultaneously at 1=12.5 seconds in the lest data of Fig. 
14. Even though the controller is sdll unslabtc, the initial 
closed-loop acceleration is much improved. 7 his is dire to 
a smooth transition from open to closed-loop. To furdicr 
improve the closed-loop bcliavior, the controller lime delay 
should be reduced. This can lie achieved by eliminating the 
full time step delay that occurs in the computations of steps 
1-5 above. Since the D matrix is null, steps 3 and 4 may 
be reordered lo compute the controller slate first and then 
the control force. This eliminates the one step delay and 
results in a stable controller as shown in Fig. 15. Note that 
it Is preferred lo reorder llic computing sequence instead 
of substituting die stale equation into the control equation 
lo eliminate lime delay as the latter results in a D matrix 
which is not null thus requiring additional compulations. 


Structure 


Switch Q 


Controltor 


figure 12. Controller closed-loop mipleittenlut ion 
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Figure 13. Closed loop responses at locations 1 and 
3 with switch A closed at 1=0 
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Time 


figure 14. Closed loop responses at lorn lions I and 
I wif h swif.cli A closed af. l=12.. r ) 



Time (second*) 

I ifinr r l. r ) Closed loop ri'sponsrs at. locations I and 
t without our slop delay 


Although the coni roller is now stable, the response in- 
die .lies aliasing because of [he relatively slow sample rale. 
When increasing the sample rale, the number of compu- 
lations required u> implement the controller is usually a 
limiting factor, lienee it is desirable to limit the size of the 
controller compulations. This is often done by controller 
stale reduction, or block diagonali/alion of the n by n stale 
matrix <t> using some real transformation. An interesting 
feature of the AVA controller is its second-order structure. 
When transformed to a first order form, the second-order 
controller equations do not represent n independent stales 
but rather n/2 in<k'|>cndcnt slates. This permits the con- 
troller rate stales to Ire expressed in terms of controller po- 
silion states. For example, a one-step implicit trapezoidal 
integration results in (he rale stales » as functions of the 
position stales d as: 

•(*)= fW*)- *(*- 1)| -«(*-!) 


where k is the lime step, h = i(t) - i(t - t). The above 
equation was used U» reduce Ihe controller stale computa- 
tions for the AVA controller which pcrmilicd a speed in- 
crease from 80 Hz. to 350 llz. for the 24 stale, 8 input 8- 
oulput controller using a VAXstation 3200 computer. Fig- 
ure 16 shows the response of the system for this update 
rate. Notice the aliasing (hat occurred while sampling at 
80 Hz is now absent and that the closed-loop performance 
is improved. 

In Fig. 16 a 32 Hz mode is being excited by the 
controller that was not noticeably excited in the excitation 
phase. This 32 Hz mode has been traced lo spatial non- 
collocation of the sensors and actuators in the original 
installation. Figure 17 shows the original installation has 
the ac tua tors (thrusters) mounted in the center or the trass 
bay faces, whereas the sensors were mounted on the trass 
hay comers. Because of ibis spatial offset, lorsion-bcnding 
coupling of the truss vihrnlinns resulted in the closer! loop 
controller exciting a torsion mode at 32 Hz. The sensor 
installation was modified as shown in Fig. 17 lo prevent 
local torsional behavior from being sensed. 

This modified installation was used in (he results or 
Fig. 18. As can be seen, the 32 Hz mode is no longer 
being driven by the controller. 

In summary, the implc menial ion of the AVA controller 
using acceleration feedback is improved by: 


(a) Not calculating controller slates during the excitation 
phase 

(b) Rcscqucncing the controller compute slcps to reduce 
delay 

(c) Updating the controller as fast as possible 

(d) Improving the spatial collocation of sensors and actua- 
tors. 


Ihe following section will describe Ihe controller per- 
formance shown in Fig. 18 in more detail. 



Time (seconds) 

Figure 1C. Cloned -loop responses at locations I and 
.3 with 350 Hz sample rate 
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IIi_Slumlaiion t ml Teat Comp»f|^ n 


In determining the success of AVA controller, each 
of the nine tnrxlcs was excited independently, and the 
controller was then activated for vibration suppression. 
"ITtc lop of Fig. 19 shows the experimental open and 
closed-loop accelerometer responses for the sixth mode 
of tlic structure, 'litis mode is controlled from location 
8. For open-loop, the structure is excited using sinusoidal 
excitation at the frequency of the mode of interest for the 
first 10 seconds. For closed loop, the structure is excited in 
the same manner, and then the controller is activated for the 
duration of each test. Shown in the bottom of Fig. 19 nrc 
the open and closed-loop lime history comjniler simulations 
performed for the same conditions. Figure 20 depicts the 
experimental and simulated open and closed-loop responses 
for the eighth mode. 

Table 3 lists the finite element and experimental fre- 
quencies, opcn-l<x>p damping, predicted closed-loop damp- 
ing, and experimental closed-loop damping. The experi- 
mental and analytical frequencies for the first nine modes 
arc fairly accurate, but at higher frequencies (not shown), 
this accuracy decreases considerably. Ihc predicted and 
actual damping values achieved when implementing this 
controller arc somewhat dilTcrcnt. It is of interest to note 
that mrxlc 2 showed no damping improvement in Ihc ex- 
periment. This is due in part to the pendulum mode accel- 
eration being nearly equal and opposite to the geometric 
change in acceleration due to gravity. Since the finite el- 
ement model of the structure d<x:s not include gravity, the 
experimentally sensed acceleration for this nuxlc is signif- 
icantly lower than predicted. Although the damping ob- 
tained with die AVA controller is not predicted well, sta- 
bility is conserved under modeling errors. 


Oiigrtial Irtr.Ult.itinn 


McxMit‘d In ‘•.fallal (on 


A< t«Ur.ilb>n 



f hrusi#t3 




Time (seconds) 

figure 18 . (’lose. | loop responses at locations I and 
willi modified srnsor insinuation 



Time (seconds) 


figure I ft. Open an. I , lose. I 
8 for mode 6 


loop responses at. location 



T ime ( seconds ) 

1 iftiirc 20. Open and closed loop responses at location 

2 for mode 8 


I i^urc 17. IVslOed .srnsor/ar tualfir installation 
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Thblc 3: Comparison of Experimental ami Analytical Frequency and Damping Results 


Mode 

Number 

I Frequency (Hz) 

F Damping 1 %) 1 

Analysis 

Experiment 

Experiment 

Open-loop 

Analysis 

Closed-loop 

Experiment 

Closed-loop 

1 

.146 

.14 

4.7 

8.7 

10.6 

2 

.148 

.14 

7.0 

10.9 

6.8 

3 

.154 

.15 

7.0 

13.0 

13.5 

4 

.720 

.73 

1.5 

1.5 

8.2 

5 

.741 

.72 

1.2 

6.5 

6.2 

6 

.868 

.90 

.6 

3.7 

4.7 

7 

1.43 

1.50 

.4 

4.6 

9.3 

8 

1.66 

1.72 

.66 

7.8 

13.1 

9 

1.83 

1.90 

.50 

4.8 

8.0 


Cvnclusiona 

The implementation of the active vibration absorber has 
been experimentally demonstrated. Using only accelera- 
tion feedback, a second order controller design is presented 
which represents an active tuned spring-mass-damper as- 
sembly, Experiment and simulation results show the con- 
troller performs well even in the presence of modeling er- 
rors. Since acceleration sensors arc relatively inexpensive 
and provide an inertial measurement, this type of sensor is 
likely to be used in structural control applications. The con- 
troller design presented herein permits direct use of accel- 
eration signals without the need for prcfiltcring. Moreover, 
the active vibration absorber has high stability robustness 
as only its performance, not its stability, is model-based. 
General implementation issues arc discussed to insure the 
render can properly use this technology in other applica- 
tions. The effects of poor implementation arc shown using 
experimental data. With proper implementation, the con- 
troller successfully augments the damping of the structure. 
This simple controller design, due to its high level of sta- 
bility robustness, has great potential for distributed, inner- 
and outer-loop control systems for spacecraft structures. 
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